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1. Introduction 



Staggered fermions had long been perceived as disadvantaged compared to Wilson fermions 
regarding the index theorem connection between (would-be) zero-modes and gauge field topology. 
For Wilson fermions, the would-be zero-modes can be identified as eigenmodes with low-lying 
real eigenvalues; these can be assigned chirality ±1 according to the sign of y'YsY' thereby de- 
termining an integer-valued index which coincides with the topological charge of the background 
lattice gauge field in accordance with the index theorem when the gauge field is not too rough 
M, ||, p. It coincides with the index obtained from the exact chiral zero-modes of the overlap Dirac 
operator [Q. In contrast, for staggered fermions, no way to identify the would-be zero-modes was 
known. They appeared to be mixed in with the other low-lying modes (all having purely imaginary 
eigenvalues) [[]], ||] and only separating out close to the continuum limit It seemed that, away 
from the continuum limit, the best one could have was a field-theoretic definition of the staggered 
fermion index [JlJ]. The latter had the disadvantages of being non-integer, requiring a renormaliza- 
tion depending on the whole ensemble of lattice gauge fields, and being significantly less capable 
than the Wilson fermion index of maintaining the index theorem in rougher backgrounds [|J]. 

Recently the consensus viewpoint described above was found to be incorrect: Staggered 
fermions do have identifiable would-be zero-modes away from the continuum limit, with identifi- 
able chiralities and integer-valued index satisfying the index theorem when the lattice gauge field 
is not too rough [0]. The would-be zero-modes, chiralities and index can be identified in a spectral 
flow approach based on a new hermitian version of the staggered Dirac operator, paralleling the 
spectral flow approach to the index for usual Wilson fermions [Q, Qj. 

Further developments along this line have led to a new version of overlap fermions built from 
staggered fermions in place of Wilson fermions [§]. The construction has the remarkable feature 
of reducing the 4 fermion flavors described by the staggered fermion to 2 flavors for the staggered 
overlap fermion. It turns out that underlying this construction is a new Wilson-type fermion, ob- 
tained by adding a Wilson-type term to the staggered fermion, which gives masses ~ 1 ja to 2 of the 
flavors while leaving the remaining 2 flavors massless. Other Wilson-type terms are also possible; 
another one which reduces the flavors from 4 to 1 was subsequently proposed in [ph. Numerical 
investigations of the 2-flavor staggered overlap fermion have been reported in [|j"0|]. The methods 
of [^, ||] were later applied to naive fermions and minimally doubled fermions [11]. 

A posteriori, these results and constructions can superficially seem quite straightforward. But 
a priori the odds were very much against any of this working out in a sensible way. There were a 
number of surprises and unexpected aspects, and these will be highlighted in the present review. 



2. Would-be zero-modes and index of the staggered Dirac operator 

In the continuum setting, the spectral flow perspective on the index of the Dirac operator D 
arises by considering the eigenvalues {A(m)} of the hermitian operator 

H{m) = Y 5 {D-m) (2.1) 

The spectral flow is defined as the net number of eigenvalues X(m) of H(m) that cross the origin, 
counted with sign ± depending on the slope of the crossing, as m is varied over some range. It can 
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be shown that the spectral flow of H{m) comes entirely from eigenvalue crossings at m = and 
equals minus the index of D. 

In the lattice setting with Wilson fermions, the spectral flow perspective |3|] is based on the 
hermitian lattice analogue of (|Q|): 

H w (m) = Y 5 (D w -m) (2.2) 

where Dw is the Wilson Dirac operator. The eigenvalue crossings of Hw(m) are in one-to-one 
correspondence with real eigenvalues of Dw, and the index of Dw (obtained from the would- 
be zero-modes, i.e. the eigenmodes with low-lying real eigenvalues) coincides with minus the 
spectral flow of the low-lying eigenvalue crossings of Hw{m). Numerical results illustrating this 
can be found, e.g., in [gj. An illustration in the d = 2 case is given in Fig.l where the eigenvalues 
of Hw (m) are plotted as functions of m. 

In the case of staggered fermions, the staggered Dirac operator D st is anti-hermitian and there- 
fore all its eigenvalues are purely imaginary. Hence the identification of would-be zero modes and 
index in the Wilson case does not carry over to the staggered case: there are no real eigenval- 



ues, and in fact the staggered analogue of (2.2), T${D st — m), is not even hermitian. The lack of 
any obvious way to distinguish the would-be zero-modes from the other low-lying eigenmodes of 
D st gave rise to the consensus viewpoint that staggered fermions are disadvantaged in this regard 
relative to Wilson fermions. 

However, it turns out that there is an alternative spectral flow approach in the staggered case 
[j7|]. Note that in the continuum setting, instead of ( |2.1| ) one can just as well use the hermitian 
operator H (m) = iD — my*, for the spectral flow perspective on the index. But now the staggered 
analogue, 

H st (m) = iD st - mF 5 (2.3) 

is also hermitian and so its spectral flow can be considered as well. Here F5 is the analogue of 
y 5 in the staggered formulation; it is hermitian and corresponds up to 0(a 2 ) discretization errors 



to 75 ® 1 in the spin(g>flavor interpretation [12]. Since H st (0) = iD st , the would-be zero-modes 
of are able to be identified as the eigenmodes with eigenvalues —iX = —iX(0) for which the 
associated flow A (m) crosses zero at a low-lying value of m. Furthermore, the sign of the slope of 
the crossing is minus the chirality of the would-be zero-mode, and hence the index is minus the 
spectral flow of H st (m) coming from the crossings at low-lying values of m. See [Q] for the details 
of this identification. 

This way of identifying the would-be zero-modes of D st from the low-lying eigenvalue cross- 
ings of H st (m) relies on an implicit assumption, namely that there is a clear separation between 
the low-lying and high-lying crossings. Actually, there is no a priori reason to believe that this 
assumption is true, even in smooth gauge field backgrounds or in the free field case. In fact one 
would expect that it is not true. The clear separation between low-lying and high-lying crossings in 
the Wilson case (as seen in Fig.l) relies crucially on the property t| = 1. But the staggered version 
F$ does not have this property. The eigenvalues of F$ are not ± 1 but are distributed throughout 
the interval [—1,1]. E.g. is an eigenvalue of F5 in the free field case; this can be seen from the 
free field momentum representation of which is riv cos2 (Pv)- In light of this one would expect 
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that, even in the free field case, the eigenvalue crossings of H st (m) will be an arbitrary mess with 
no clear separation into low-lying and high-lying crossings. 

The first and biggest surprise in all this - a miracle in fact - is that, contrary to expectations, the 
spectral flow of H st (m) does have a clear separation between low-lying and high-lying eigenvalue 
crossings, at least when the gauge field is not too rough. In fact in the free field case there are no 



high-lying crossings at all, cf. the bound ( [2.7| ) below. Fig. 2 shows the spectral flow in a moderately 
roughened U(l) background with topological charge Q = 1 on a 2-dimensional lattice. Now there 
are high-lying crossings, but they are clearly separated from the low-lying ones, so the would-be 
zero-modes of D st , their chiralities, and index, can be unambiguously identified. 

The absence of high-lying eigenvalue crossings for H st (m) in the free field case can be seen 
analytically as follows. A simple calculation of H st (m) 2 in the free field momentum representation 
gives 

HstH 2 = £sin 2 (/^) + m 2 ncos 2 (/? v ) (2.4) 
H v 

Set = sm(pn), c v = cos(/? v ). Then, in the case of 2 spacetime dimensions, starting from 
H st (m) 2 = s\ + sl + m 2 {\ - s\)(\ -s$), we find 

H st (m) 2 =m 2 + (l-m 2 )(s 2 +sl)+m 2 s 2 l s 2 , > m 2 for < \m\ < 1 (2.5) 

and 

H st (m) 2 = l+s 2 l s 2 2 + (m 2 -l){\-s\){\-s 2 2 ) >1 for |m| > 1 (2.6) 

Note that both of these bounds are saturated. Identical bounds can be derived in the d = 4 case, 
although the derivations are more complicated. Hence in the free field case, for both d = 2 and 
d = 4 dimensions (and probably also for higher dimensions), we have 

. , 9 I m 2 for \m\ < 1 

HAmffree > { T f 7 , (2-7) 

J I 1 for \m\ > 1 

This bound has a generalization to the case of gauge fields satisfying an "admissibility condition" 
on the plaquettes that implies a separation between low-lying and high-lying eigenvalue crossings 
for H st (m) when the e in the condition is sufficiently small - see [Q] for details. Thus the situation is 
analogous to the Wilson case where the admissibility condition guarantees the separation between 
low-lying and high-lying crossings for the hermitian Wilson operator Hw(m) fll3[]. 

Comparing Fig.'s 1 and 2 we see that the form of the spectral flow of H st (m) is very different 
from the Wilson case, and the separation between the low-lying and high-lying crossings is much 
larger. However, this is not one of the surprises alluded to in the abstract. Instead, the surprise here 
is that there is no surprise - the staggered spectral flow has the same form as in the Wilson case 
once the correct interpretation of the hermitian staggered operator H st (m) is identified. It turns out 
that the parameter m in the staggered case should be identified not with the corresponding m in 
H\v(m) but with the Wilson parameter r in the Wilson case. This will be explained further below; 
see Fig. 3. 
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Since the staggered fermion in d dimensions describes 2 d l 2 flavors, the index theorem in this 
case should be index(D it ) = 2 d / 2 (—\) d / 2 Q. This is confirmed by numerical results in smooth 
enough backgrounds. E.g. in Fig.2 the two positive slope low-lying crossings in the Q = 1 back- 
ground in 2 dimensions imply that the index is —2 in accordance with the index theorem. The 
eigenvalues of D st in this background correspond to m = in Fig.2. The would-be zero-modes of 
D st can be identified as the two eigenmodes with the low-lying eigenvalues belonging to the two 
eigenvalue flows that cross the origin. 




-8 -4 -2 2 4 6 -10 -5 5 10 



Figure 1 : Eigenvalue flow of Hw (m) in a Figure 2: Eigenvalue flow of H st (m) in the 

Q=l background on a 2-dim lattice. same Q=\ background as Fig. 1 . 



3. Staggered overlap construction 

In the Wilson case, the spectral flow perspective on the index leads to 

index(D w ) = _i Tr Hw ^ m \ = index(D ov ) (3.1) 



for any mo in the region between where the low-lying and high-lying eigenvalue crossings of 
H w (m) occur (e.g. m = 1/a) [pi Bfl. Here D ov = -(l + y 5 " wim) is the overlap Dirac op- 

" V y/H w (m ) 2 J 

erator The intimate connection between the Wilson index, the hermitian operator Hy/{m) and 
the overlap Dirac operator suggests there may exist a staggered version of the overlap Dirac op- 
erator connected to the staggered index and staggered hermitian operator H st (m) discussed above. 
But there is a problem: The connections and properties of the overlap Dirac operator in the Wilson 
case rely crucially on the property y| = 1. (E.g. without it the GW relation would not hold and 
D ov would not have exact zero-modes.) Therefore the natural replacement 75 — > F$ is not possible 
when constructing the overlap operator in the staggered case since r| 7^ 1. Attempts to construct 
different versions of 75 in the staggered case which do satisfy y 2 = 1 invariably lead to unnatural, 



problematic operators which violate either lattice rotation invariance or gauge invariance QL4J. 

This 'Tj 7^ 1 problem" initially appears insurmountable, but the second main surprise in all 
this is that it does have a solution The theoretical idea behind the solution is as follows. In 
the staggered setting there is a naturally arising operator which squares to the identity, namely 
T55, acting on the staggered fermion fields by Y^x{ x ) = (— l)" 1+ "' + ' v z( x )- H nas the spinigiflavor 
interpretation 75 (g> 75, which is not what we want. But if the staggered overlap construction can be 
set up such that the physical flavors are those with positive flavor-chirality under 1 ® 75 then 75 <g) 75 
will be the same as 75 ® 1 on the physical flavors, and then 1^55 may be used for the role of 75 in the 
staggered overlap construction. 
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In fact this can be achieved simply by replacing 75 — > Y55 and Hw — >■ H st in the overlap formula 
for D ov . The key observation is that in this construction we have Y^Hst^mo) = iF^D^ — moY^^Y^ 
and T55r5 has the spinigiflavor interpretation (75 (g) 75) (75 <8> 1) = 1 ® 75 up to 0(a 2 ) effects. Thus 
the 2 fermion flavors with positive flavor-chirality get a negative mass from the — 7^(^55^ term, 
and hence become massless flavors of the staggered overlap fermion, while the 2 flavors with 
negative flavor-chirality get positive masses from this term and hence become heavy, decoupling 
flavors of the staggered overlap fermion with masses ~ I /a, just like the "doubler" species in the 
usual overlap fermion construction. Note that the exact flavored chiral symmetry {D st ,Y^} = 
of the staggered fermion hereby becomes an exact unflavored GW chiral symmetry {D sov , ^5 } = 
aD sov rs5D sov of the resulting staggered overlap Dirac operator D sov . Moreover, a staggered version 
of the index relations ([U]) holds [0, §: ±index(D,,) = -^Tr-^pL = index(D, 0V ). The factor \ 

multiplying index(D 5 ,) reflects the reduction from 4 to 2 flavors in the staggered overlap fermion. 

The interpretation of the staggered overlap fermion becomes more straightforward if we change 
the hermitian staggered operator by H st {m) — > Y=,sD st — mY^ = Y^{D st — mY^Ys). As mentioned 
in []|], this operator is closely related to, and has the same eigenvalue spectrum as, the previous 
operator iD st —mY^. Everything in the preceding continues to hold with this new H st (m). The 
staggered overlap Dirac operator takes a more recognizable form though: it can now be written as 

D sov = - (l + (D st - mo Y 55 Y 5 )^= ===== = = = ) (3.2) 

From this we see that underlying the staggered overlap construction is a new staggered version of 
Wilson fermions with the Dirac operator 

D sW =D st + W st , W« = -(l-r 55 r 5 ). (3.3) 

a 

The "Wilson term" W st decouples the negative flavor-chirality modes by giving them mass 2r/a 
while keeping the two positive flavor-chirality modes as the physical modes. Hence D s w describes 
two physical quark flavors on which Ts = ^5 up to 0(a) effects. It has the hermiticity d] w = 
^5sD s wTs5. A 2-flavor overlap fermion can then be obtained by taking D s w —m with m = 
p G (0, 2) as the kernel in the usual overlap construction. For ^ = niQ and p = 1 this is precisely 



the 2-flavor staggered overlap Dirac operator D sov obtained above in (3.2). But now we see that it 
can be generalized to any p G (0,2). Furthermore, the role of the parameter mo in the staggered 
overlap construction is hereby clarified: it is analogous to the Wilson parameter in the usual overlap 
construction. 



The general staggered overlap operator can also be expressed as D sov = -(1 +1^55 



where H s w is another hermitian staggered operator given by 1 

H sW (m) = Y 55 (D sW -m) = Y 55 D st - Y 5 + (1 - m)Y 55 (3.4) 

This is the true analogue of the hermitian Wilson operator Hyy(m) = 75 (£V — m )> an d its spectral 
flow has a similar form (Fig.4). 2 On the other hand, the spectral flow of the hermitian Wilson 

'We set r = 1 and use lattice units to get the second equality. 

2 In 2 dimensions the staggered Wilson fermion has one physical flavor and one doubler whereas the usual Wilson 
fermion has three doublet's. This explains why there is one high-lying eigenvalue crossing in Fig.4 and three high-lying 
crossings in Fig. 1 . 
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operator as a function of the Wilson parameter r, with fixed m = 1 (Fig.3), has a similar form to 
the spectral flow of our previous hermitian staggered operator H st (m) as a function of m (Fig.2) as 
anticipated. 3 












Figure 3: Eigenvalue flow of Hw(m = I) 
as a function of the Wilson parameter r. 



Figure 4: Eigenvalue flow of H s w(m) in 
the same Q—l background as Fig.s 1,2,3. 



To summarize, three new lattice fermion formulations have been introduced: staggered ver- 
sions of Wilson fermions, domain wall fermions 4 and overlap fermions. The relation of ordinary 
staggered fermions to these staggered-based formulations is analogous to the relation of naive 
fermions to the Wilson-based formulations. The crucial property of the "Wilson term" in the stag- 
gered version of Wilson fermions is that it splits the 4 flavors into pairs with positive and negative 
flavor-chirality, giving a mass ~ I /a to the latter while leaving the former pair as the massless phys- 
ical fermions. This allows ^5 ~ 75 <g> 75 to be used in place of Y=,, thus overcoming the 'Tg 7^ 1 
problem". The staggered versions of domain wall and overlap fermions can then be constructed 
simply by replacing 75 — > r 5 $ and D w — > D sW in the usual formulations. 

Each of the staggered-based formulations is a new alternative and competitor to the corre- 
sponding Wilson-based formulation, and should be more computationally efficient since the con- 
structions start from staggered rather than naive fermions. However, the gain in efficiency is re- 
duced because the staggered "Wilson term" is less local; it involves Ts which is a 4-link operator 



[10]. This reduction in efficiency may possibly be ameliorated by smearing the links 

Much remains to be done to clarify the theoretical and practical properties of the new staggered- 
based formulations. They break some of the symmetries of the original staggered fermion and the 
consequences of this need to be investigated and clarified. This may be done in lattice perturba- 
tion theoiy and also non-perturbatively by numerical lattice QCD calculations. E.g. it should be 
checked for the staggered Wilson fermion that a chiral (massless) limit can be reached by tuning 
the bare quark mass. 5 For the staggered domain wall fermion and staggered overlap fermion it 
should be checked that the lattice QCD theory has an approximately and exactly massless phase, 
respectively. Work on this is currently underway. 
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